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Abstract 

Based on an isospin- and momentum-dependent hadronic transport model IBUU04, effects of the 
nuclear symmetry energy on the single and double 7r“/7r+ ratios in central reactions of ^^^Sn-|-^^^Sn 
and ^^^Sn-|-^^^Sn at a beam energy of 400 MeV/nucleon are studied. It is found that around 
the Coulomb peak of the single vr^/vr"’' ratio the double ratio taken from the two isotopic 

reactions retains about the same sensitivity to the density dependence of nuclear symmetry energy. 
Because the double ' k ~ ratio can reduce significantly the systematic errors, it is thus a more 
effective probe for the high-density behavior of the nuclear symmetry energy. 
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I. INTRODUCTION 


After about three decades of intensive efforts in both nuclear experiments and theories, the 
equation of state (EOS) of isospin symmetric nuclear matter is now relatively well determined 


mainly by studying collective flows in heavy-ion collisions ^ and nuclear giant monopole 
resonances The major remaining uncertainty about the EOS of symmetric nuclear 
matter is due to our poor knowledge about the density dependence of the nuclear symmetry 


energy [1|,|3|,|^. Therefore, the new challenge is to determine the EOS of isospin asymmetric 
nuclear matter, especially the density dependence of the nuclear symmetry energy. Besides 
the great interests in nuclear physics, the EOS of asymmetric nuclear matter is also crucial 
in many astrophysical processes, especially in connection with the structure of neutron stars 


and the dynamical evolution of proto-neutron stars p. Fortunately, heavy-ion reactions, 
especially those induced by radioactive beams, provide a unique opportunity to constrain 


the EOS of asymmetric nuclear matter 


In fact, considerable progress has been made 


recently in determining the density dependence of the nuclear symmetry energy around the 
normal nuclear matter densit y fr om studying the isospin diffusion in heavy-ion reactions at 


intermediate energies 


HQ. 


However, much more work is still needed to probe the 
high-density behavior of the nuclear symmetry energy. 

A crucial task is to hnd experimental observables that are sensitive to the density depen¬ 
dence of the nuclear symmetry energy. A number of such observables have been already iden- 
tihed in heavy-ion collisions induced by neutron-rich nuclei, such as the free neutron/proton 
121, the isospin fractionation 13|, [ij, [l^, |ly, [iT], ll^ , the neutr on-p roton transverse 


ratio 


differential flow 


the neutron-proton correlation function 


23, t/3He |23l |2|, 


the isospin c 


29 


30, 


iffusion j25|, 123, the proton differential elliptic flow 22| and the vr ratio 


31|. Generally, the long range Coulomb force on charged particles plays an 


important role in the above observables. It is thus important to distinguish the effects due 
to the symmetry potentials from those due to the Coulomb potentials. It is also useful 
to understand the interplay between these two kinds of potentials. Moreover, to extract 
accurately information about the symmetry energy one has to reduce as much as possible 
the systematic errors involved in the observables used in experiments. For this purpose one 
normally studies the ratios or relative values of experimental observables from two reaction 
systems using different isotopes of the same element. The hrst theoretical investigation of 
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such an observable, the double neutron/proton ratio of pre-equilibrium nucleons, has been 
made recently in Ref. Since the single ti~ ratio in heavy-ion collisions induced by 


neutron-rich nuclei has been s’ 
nuclear symmetry energy |^. 


lown to be a useful probe of the high-density behavior of the 


30 


3 , 


we study here the effects of the symmetry energy on 


the double tt /tt’*' ratio from the reactions of and ^^^Sn-|-^^^Sn, i.e., the ratio 

of 7r“/7r+ from ^^^Sn-|-^^'^Sn over that from ^^^Sn-|-^^^Sn, using the IBUU04 model. It is 


well known that the single vr ratio has a Coulomb peak at certain pion 


inetic energy 


361. It is thus 


depending on the system and the impact parameter of the reaction 33|, |3^ |35 
especially interesting to examine the sensitivity of the double 'k ~ratio to the symmetry 
energy around the Coulomb peak. We hnd that around the Coulomb peak the double 'k~ 
ratio has about the same sensitivity to the symmetry energy as the single ti~ ratio while 
having the advantage of reduced systematic errors. 

II. A BRIEF INTRODUCTION TO THE IBUU04 TRANSPORT MODEL 


In the IBUU04 model, besides nucleons, A and iV* resonances as well as pions and 
their isospin-dependent dynamics are included. The initial neutron and proton density 
distributions of the projectile and target are obtained by using the relativistic mean field 
theory. The experimental free-space nucleon-nucleon (NN) scattering cross sections and 
the in-medium NN cross sections can be used optionally. In the present work, we use the 
isospin-dependent in-medium NN elastic cross sections from the scaling model according 
to nucleon effective masses ^|. For the inelastic cross sections we use the experimental 
data from free space NN collisions since the in-medium inelastic NN cross sections are still 
very much controversial. The total and differential cross sections for all other particles are 
taken either from experimental data or obtained by using the detailed balance formula. The 
isospin dependent phase-space distribution functions of the particles involved are solved by 
using the test-particle method numerically. The isospin-dependence of Pauli blockings for 
fermions is also considered. More details can be found in Refs. |ll[ 1^1^. 1^ 1^ . The 
momentum-dependent single nucleon potential (MDI) adopted here is 0 


f/(p,5,p,r) = A^{x)—+ Ai{x) — 

Po Po 
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In the above equation, 6 = (Pn—Pp) / (Pn+Pp) is the isospin asymmetry parameter, p = Pn+Pp 
is the baryon density and Pn,Pp are the neutron and proton densities, respectively, r = 
l/2(—1/2) for neutron (proton) and t ^ t', a = 4/3, /r(r, p) is the phase-space distribution 
function at coordinate r and momentum p. The parameters A„(a;), A;(a;), 5, a,r, a,r' 
and A were set by reproducing the momentum-dependent potential U (p, 6, p, r) predicted 
by the Gogny Hartree-Fock and/or the Brueckner-Hartree-Fock calculations, the saturation 
properties of symmetric nuclear matter and the symmetry energy of about 32 MeV at normal 
nuclear matter density po = 0.16 fm“^. The incompressibility of symmetric nuclear matter 
at normal density is set to be 211 MeV. The parameters Au{x) and Ai{x) depend on the x 
parameter according to 


Au{x) = -95.98 


25 25 

-X, Ai{x) = —120.57 H- rx, 


( 2 ) 


a + 1 a + 1 

where B = 106.35 MeV. A = is the nucleon Fermi momentum in symmetric nuclear 
matter. Cry = —103.4 MeV and a,r = —11.7 MeV. The Cry and a,T terms are the 
momentum-dependent interactions of a nucleon with unlike and like nucleons in the sur¬ 
rounding nuclear matter. The parameter x is introduced to mimic various density-dependent 
symmetry energies Esym{p) predicted by microscopic and phenomenological many-body ap¬ 


proaches H- The isoscalar part {Un + Up)/2 j37|, |39| of the single nucleon potential was 


shown to be in good agreement with that of the variational man y-b ody calculations 


m 


and the results of the BHF approach including three-body forces 
{Un — Up)/25 37| is consistent with the experimental Lane potential 
According to essentially all microscopic model calculations, see e.g., j4 
isospin asymmetric nuclear matter can be expressed as 


. The isovector part 


the EOS for 


E{p, 5) = E{p, 0) + aym(p)^' + 0(h"), 


(3) 


where E(p, 0) is the energy per nucleon of symmetric nuclear matter, and aym(p) is the 
nuclear symmetry energy. With the single particle potential [/(p, 5, p, r), for a given value 
X, one can readily calculate the symmetry energy Egy^iip) as a function of density. Noticing 
that the isospin diffusion data from NSCL/MSU have constrained the value of x to be 
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FIG. 1: (Color online) Density dependence of nuclear symmetry energy using the MDI interaction 


with X = 0 and x = — 1. Taken from 




between 0 and —1 for nuclear matter densities less than about 1.2po OL M 1j|. in the present 
work we thus consider only the two values of x = 0 and x = — 1. Shown in Fig. ^ is the 
density dependence of the nuclear symmetry energy with the two x parameters. It is seen 
that the case of x = 0 gives a softer symmetry energy than that of x = —1 and the difference 
becomes larger at higher densities. 


III. RESULTS AND DISCUSSIONS 


It was found earlier that the single tt ~ratio from heavy-ion collisions induced by 
neutron-rich nuclei can be used to probe the high density behavior of the nuclear symmetry 
energy j3, Q 1^ In Fig. O we show the kinetic energy distribution of the single 
'n~ j-K^ ratio for the reactions of and ^^^Sn-|-^^^Sn at a beam energy of 400 

MeV/nucleon and an impact parameter of 6 = 1 fm with the stiff (x = —1) and soft 
(x = 0) symmetry energy, respectively. In order to obtain good statistics, we used 12000 
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events for each reaction in the present work. It is seen that the overall magnitude of 'k~ 
ratio is larger for the neutron-rich system ^^^Sn-|-^^'^Sn than for the neutron-dehcient system 
^^^Sn-|-^^^Sn as expected. For the reaction ^^^Sn-|-^^^Sn the single tt~/tt^ ratio is not so 
sensitive to the symmetry energy due to the small isospin asymmetry. However, for the 
neutron-rich system ^^^Sn-|-^^‘^Sn the single ratio is sensitive to the symmetry energy 


and this is consistent with previous studies 2^ |29|, |30 , . It is further seen that the soft 

symmetry energy {x = 0) leads to a larger single ti~ ratio than the stiff one {x = —1). 
This is mainly because the high density region where most pions are produced are more 
neutron-rich with the softer symmetry energy as a result of the isospin fractionation 0 



FIG. 2: (Color online) Kinetic energy distribution of the single vr^/vr'*' ratio for ^^^Sn-|-^^^Sn and 
^^^Sn-|-^^^Sn at a beam energy of 400 MeV/nucleon and an impact parameter of 6 = 1 fm with the 
stiff {x = —1) and soft (x = 0) symmetry energies. The dash-dotted line is the single Tr^/Tr"*" ratio 
obtained by turning off the Coulomb potentials in the ^^^Sn-|-^^^Sn reaction. 

From Fig. |21 it is also interesting to see that the single 7r~/7r~^ ratio exhibits a peak at a 
pion kinetic energy of about 45 MeV in all cases considered here. In order to understand the 
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origin of this peak, we also calculated the single 'k ~ratios in both reactions by turning 
off the Coulomb potentials for all charged particles. As an example, shown in Fig. |21 with the 
dash-dotted line is the single tt~ ratio obtained by turning off the Coulomb potentials in 
the ^^^Sn-|-^^^Sn reaction. It is seen that the single ratio now becomes approxim ately 

a constant of about 2.4. The latter is what one expects based on the A resonance model (46j |. 
According to the latter the 7r“/7r+ ratio is approximately + NZ)/+ NZ) Ri {N/ZY 
in central heavy-ion reactions with N and Z being the total neutron and proton numbers in 
the participant region |^. For central ^^^Sn-|-^^^Sn reactions, the value is about 2.43. At 
400 MeV/nucleon, pions are almost exclusively produced via the A resonances j^, it thus 
should not be a surprise to see the agreement with the A resonance model expectation. The 
comparison of calculations with and without the Coulomb potentials indicates clearly that 
the peak observed in the single 7r“/7r+ ratio is indeed due to the Coulomb effects. The ti~ 
ratio carries some information about the symmetry energy mainly because it is sensitive to 
the isospin asymmetry of the nucleonic matter where pions are produced. This information 
might be distorted but is not completely lost because of the Coulomb interactions of pions 
with other particles. It is thus natural to look for signals of the symmetry energy in kinematic 
regions where the vr^/vr^ ratio reaches its maximum. In this regard, the Coulomb peak is 
actually very useful for studying the symmetry energy. It is necessary to stress that in many 
situations the Coulomb peak will simply appear at zero instead of a finite kinetic energy. 
One would then need to concentrate on the tt ~ratio of low energy pions. In addition, it 
should be mentioned that most of pions are produced in the high density nucleonic matter 
(about 2po) through the A resonances and thus carry important information about the high 
density behavior of the symmetry energy. Since the pions at lower kinetic energies around 
the Coulomb peak experience many rescatterings with nucleons at high or low densities and 
the charged pions considered here also feel the Coulomb potential from protons at different 
densities, the information on high density symmetry energy from the lower en ergy pions may 
be distorted partially by the low density behavior of the symmetry energy [2l|. However, 
in the present work, the soft {x = 0) and stiff (x = —1) symmetry energies have slight 
difference at low densities and the large difference appears at high densities (about 2po) as 
shown in Fig. ^ and we thus expect the observed symmetry energy effects on the energy 
dependence of the ratio mainly reflect (though not completely) information on the 

high density behavior of the symmetry energy. 
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FIG. 3: (Color online) Kinetic energy dependence of the double 'k~ j'K^ ratio of ^^^Sn+^^'^Sn over 
^^^Sn+^^^Sn at a beam energy of 400 MeV/nucleon and an impact parameter 6 = 1 fm with the 
stiff (x = — 1) and soft (x = 0) symmetry energies. 

To reduce the systematic errors, it is more useful to study the double 'k ~ratio in 
the reactions of and A practically critical question is whether the 

sensitivity to the symmetry energy observed in the single ti ~ratio can be sustained by 
the double ratio. To answer this question we examine in Fig. El the double 'k ~ratio for 
the two reactions considered. It is seen that the kinetic energy dependence of the double 
'K ~ratio is rather different for the stiff (x = — 1) and soft (x = 0) symmetry energies. 
The double 7r“/7r+ ratio is quite flat for x = 0 while displaying a concave structure for 
X = — 1 around the Coulomb peak. These different behaviors can be understood from the 
corresponding single tt~ ratios in the two reactions as shown in Fig. |21 It is reassuring to 
see that around the Coulomb peak the double tt~/ti^ ratio is still sensitive to the symmetry 
energy. Compared with the single tt~ ratio, the kinetic energy dependence of the double 
'K ~ratio becomes weaker. This is because the effects of the Coulomb potentials are 
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reduced in the double 7r“/7r+ ratio. We note that the double 7r“/7r+ ratio displays an 
opposite symmetry energy dependence compared with the double n/p ratio for free nucleons 
shown in Ref. jS^. This is understandable since the soft symmetry energy leads to a 
neutron-richer dense matter in the heayy-ion collisions induced by neutron-rich nuclei and 
thus more 7r“’s would be produced due to more neutron-neutron inelastic scatterings. Since 
the soft symmetry energy leads to a neutron-richer dense nucleonic matter, the n/p ratio 
for free nucleons is therefore expected to be smaller due to the charge conseryation. 


IV. SUMMARY 

In summary, using the isospin- and momentum-dependent hadronic transport model 
IBUU04, we studied the single and double /ti^ ratios and their dependence on the nuclear 
symmetry energy in the central reactions of ^^^Sn-|-^^^Sn and ^^^Sn-|-^^^Sn at a beam energy 
of 400 MeV/nucleon. We found that the double 'k ~ratio retains the same sensitiyity to 
the symmetry energy as shown by the single tt~/ti^ ratio around its Coulomb peak in the 
neutron-richer system inyolyed. These features are useful for extracting information about 
the nuclear symmetry energy from future experimental data. Because the double 'n~ 
ratio can reduce signihcantly the systematic errors, it is thus a more effectiye probe for the 
high-density behayior of the nuclear symmetry energy. 
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